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^"""^ , Quantum shape-phase transitions in odd-even nuclei are investigated in the 

1 framework of the interacting boson-fermion model. Classical and quantum anal- 

ysis show that the presence of the odd fermion strongly influences the location 
and nature of the phase transition, especially near the critical point. Experi- 
mental evidence for the occurrence of spherical to axially-deformed transitions 
in odd-proton nuclei Pm, Eu and Tb (Z=61, 63, 65) is presented. 
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■ Quantum phase transitions (QPTs) are qualitative changes in the 

ground state properties of a physical system induced by a variation of pa- 
rameters in the quantum Hamiltonian. These structural modifications have 
found a variety of applications in diverse areas of physics. 1 One of these 
applications is to atomic nuclei, where QPTs in even-even nuclei have been 
extensively investigated (for a review, see 2-4 ) within the framework of the 
Interacting Boson Model (IBM), 5 a model of nuclei in terms of correlated 
pairs of valence nucleons with L = 0, 2 treated as (s, d) bosons. 5 In the 
present contribution we extend these studies to odd-even nuclei making use 
of the Interacting Boson Fermion Model (IBFM), 6 incorporating an addi- 
tional unpaired fermion with angular momentum j. The reported results 
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portray general properties of QPTs in mixed Bose-Fermi systems 7 and are 
illustrated for j = 11/2. 

We consider the Hamiltonian of a system of N monopole (s T ) and 
quadrupole (djj bosons and a single- j fermion (a] m ) 

H = H B + Hf + Vbf , (1) 



with 





= £a 




(2a) 


H F 


— Ej 71 j , 


(2b) 


Vbf 


= rg x -(a]a,) (2) +A v /2j + l 


: [(da j )W(da))W]<- '> : . (2c) 


Here ha 


= <ft 


• d and fij are the (i-boson and fermion number operators 



respectively and Q x = (d^s + s^d)^ + x(d^d)( 2 \ The parameter e is the 
scale of the boson energy, Sj is the energy of the single fermion, T and A 
are, respectively, the strengths of the quadrupole and exchange Bose-Fermi 
interactions. QPTs of the purely bosonic part of the Hamiltonian Hb have 
been extensively investigated. 2-4 There are two control parameters £ and 
X- For fixed \-> as one varies £, < £ < 1, the bosonic system undergoes a 
QPT. The phase transition is first order for \ 7^ and becomes second order 
at x — 0. No phase transition occurs as a function of x- In this contribution, 
we take x = m which case the spherical phase has U(5) symmetry 

(£ = 0) and the axially-deformed phase has SU(3) symmetry (£ = l). 5 The 
critical point, separating the two phases occurs at £ c = 1/2. 

A complete study of the properties of quantum phase transitions neces- 
sitates both a classical and a quantal analysis. A classical analysis amounts 
to constructing the combined Bose-Fermi potential energy surface (Landau 
potential) and minimizing it with respect to the classical variables. To this 
end, we introduce a boson condensate 5 

|JV;/3, 7 ) = (M)- 1/2 K(A7)] JV |0) , (3) 

b\ = (1 + ,3 2 )- 1/2 [,3 cos 7 (4 + /3 sin 7(4 + 4 2 )/\/2 + s% in terms of the 
classical variables f3, 7. By integrating out the boson degrees of freedom, 
i.e., by taking the expectation value of H (1) in the boson condensate, 
H(N;f3,7) — (N; f3, 7|.ff|iV; /?, 7), one obtains the fermion Hamiltonian 

9mi,rri2 

(AT;/3,7)4 roi a jiro2 , (4) 

mi ,rri2 

where Eb = (N; f3,f\H b\N; /3, 7). The matrix g mi ,m 2 {N; (3, 7) is a real, 
symmetric matrix, which depends on the Bose-Fermi couplings, T and A. 
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Its numerical diagonalization yields the single particle eigenvalues ei((3,j) 
and eigenfunctions ipi(/3, 7), i = 1, 2, + |. These are the single-particle 
levels in the deformed ((3, 7) field generated by the bosons. For 7 = 0° the 
eigenvalues are given in explicit analytic form 8 

X K (P) = -NT I (77^2) ^ (2 - f*X\[f} PA*K 2 j(j + 1)] 

-NA^J^y2j + l)P?l3K 2 -j(3 + l)f} , (5) 

where Pj = [(2j - l)j(2j + l)(j + l)(2j + 3)] 7 . The eigenvalues can be 
labelled by the angular momentum projection on the intrinsic axis 3, K3 = 
K = and they are doubly degenerate. Similarly, the eigenvalues 

for 7 = 60° (oblate axial symmetry), are given by ujk 2 {P) = Xk 3 ^-k 2 (~ 
where K 2 is the corresponding projection on the axis 2. Once the eigenvalues 
have been obtained, one can calculate the total energy functional (Landau 
potential for the combined Bose-Fermi system) 

Ei(N; p, 7; £, x; T, A) = E B (N; /3, 7; X ) + e, + e^N; 0, 7; X ; T, A) . (6) 

Minimization of Ei with respect to /3 and 7 gives the equilibrium values 
/3 e ,7 e (the classical order parameters) for each state. In general, Eb oc iV 2 
and ei oc N, however, as shown below, contrary to naive expectations, the 
odd fermion greatly influences the phase transition. 

To facilitate the study of QPTs in the IBFM we set 9 r = — 2e ^- Qjj 
with Qjj — (j\\y^\\j) and examine the change in (/? e ,7 e ) upon variation 
of the control parameter £. For pure quadrupole coupling (A = 0), one can 
see from the left panel of Fig. 1 that the phase transition is washed out 
for states 1,2,3 and enhanced for states 4,5,6. The critical point is approx- 
imately at the same location as for the purely bosonic case (£g ~0.51) for 
state 4 (£| <~0.50), but it is moved to larger values for state 5 (£f ~0.53) 
and state 6 (£g ~0.58). The values of j e below and around the critical point 
are no longer zero. After the critical point, all states become those of a sin- 
gle particle in an axially deformed field (j e — 0°) with equal deformation 
Pe,i — Pe,B {Pe,B being the equilibrium deformation of Eb)- States in this 
region can be labelled by K — ^, |, ...-j, corresponding to states i = 1, 6. 
The equilibrium values when A = —0.287 are shown in the right panel of 
Fig. 1. We see that the modifications induced by the presence of a fermion 
to the phase transition are more dramatic when the exchange interaction is 
added, especially in the location of the critical point £f . The state 6 moves 
considerably to the right of £3. In summary, QPTs in the presence of an 
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Fig. 1. Equilibrium values, /3 e i (top part) and (bottom part) as a function of the 
control parameter £ in the U(5)-SU(3) transition (x = — ^-), and N = 10. Left: A = 0. 
Right: A = —0.287. States are labelled by the index i = 1, 6. 

odd fcrmion have a different behavior as a function of the control parameter 
£ than in the purely bosonic system. The odd fermion acts as a catalyst for 
some states and as a retardative for others. 

The quantal analysis of QPTs is done by diagonalizing numerically the 
Hamiltonian (1) for finite N. In the semi-microscopic version of the IBFM, 
the Bose-Fermi interactions are given in terms of the BCS occupation prob- 
abilities, Uj and Vj, with + Vj = 1. Specifically, T — T s 2 (Uj — u|) Qjj, 
A = — A s 8>/EujVjQ?j/(2j + 1) and we set T s = ^ e . The correlation di- 
agrams for pure quadrupole coupling (A = 0) are presented in Fig. 2. They 
describe how the energy levels evolve from the U(5) phase to the SU(3) 
phase. They also display particle-hole conjugation, that is the transforma- 
tion (nf>u). As quantal order parameters we consider here the expectation 
value of fid in the states i = 1, 6, 

(1) = (gj \n d \ gi) m 
i N K ) 

This is shown in Fig. 3 top part. Since the order parameters v\ ^ are related 
to the square of the classical order parameters (3 e ^, this figure is related to 
Fig. 1 to which it corresponds in the limit N — > oo. The derivative of 

v\ in the ground state, ^ is also shown in Fig. 3 bottom part. (This 
quantity diverges when N — > oo). From this figure one sees clearly that the 
transition is made sharper by the presence of the fermion for some states, 
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Fig. 2. Correlation diagram for a j = 11/2 Fig. 3. The quantal order parameters, 

particle coupled to a system of (s, d) bosons v\ . , Eq. (7) (top part), and its deriva- 

undergoing a U(5)-SU(3) phase transition, tive, (bottom part), for yrast states with 

Top part v 2 = 1, bottom part v 2 = 0. The J = .., 1 as a function of the control 

interaction is purely quadrupole. parameter, f . Here N = 10 and v 2 = 1. 

11/2,9/2,7/2, while is made smoother for others, 5/2,3/2,1/2, a result 
already seen in the classical analysis. 

One of the experimental signatures of QPT in nuclei is the two-neutron 
separation energies, Si n = — [Eo(N + 1) — E (N)], which can be related to 
the derivative of the ground state energy, Eq, with respect to the control 
parameter, S^n can be written as a smooth contribution linear in the 
boson number N, plus the contribution of the deformation 5 

S 2n = ~A 2n - B 2n N + S(2n) dci . (8) 

In order to emphasize the occurrence of the phase transition it is convenient 
to plot the deformation contribution only, obtained from the data by sub- 
tracting the linear dependence, as a function of N. In previous studies of 
the purely bosonic part it has been shown that N is approximately propor- 
tional to the control parameter £. 2 The experimental values of S , (2n) < j e f are 



October 5, 2011 0:43 



WSPC - Proceedings Trim Size: 9in x 6in cgsl4lev 



6 



Fig. 4. The contribution of deforma- 
tion to the two-neutron separation ener- 
gies, S(2n) < j e f for even-even 6oNd-62Sm- 
64 Gd nuclei (top) and odd-even eiPm- 
63Eu-fj5Tb nuclei (bottom), plotted as 
a function of neutron number. The con- 
tribution is enhanced in odd-even nuclei 
by approximately 300 keV (at neutron 
number 92). Also the rise between neu- 
tron numbers 88 and 90 is sharper in 
odd-even nuclei than in even-even nu- 
clei. In the limit TV — > oo (no finite size 
scaling) the quantity S(2n) ( j c f should 
be zero before the critical value and fi- 
nite and large after that. The expected 
behaviour of for the U(5)-SU(3) 

transition and N = 10 is shown in the 
inset. 

shown in the top part of Fig. 4 for even-even nuclei (purely bosonic) and 
in the bottom part for odd-even nuclei (bosonic plus one fermion). They 
are obtained from the empirical data with A 2n — —14.61, —15.82, —16.997 
MeV for Nd-Sm-Gd, respectively, and B 2n = 0.657 MeV, and with A 2n = 
-15.185, -16.37, -17.672 MeV for Pm-Eu-Tb, and B 2n = 0.670 MeV. Pre- 
cursors of the phase transition are visible in all six nuclei between neutron 
numbers 88 and 90 in both, and, most importantly, appears to be enhanced 
in odd-even nuclei relative to the even-even case. 

This work was supported in part by the U.S. -Israel Binational Science 
Foundation and in part by DOE Grant No. DE-FG-02-91ER40608. 
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